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Abstract. In this paper we consider the Brownian motion with jump boundary and present 
a new proof of a recent result of Li, Leung and Rakesh concerning the exact convergence rate 
in the one-dimensional case. Our methods are different and mainly probabilistic relying on 
coupling methods adapted to the special situation under investigation. Moreover we answer a 
question raised by Ben-Ari and Pinsky concerning the dependence of the spectral gap on the 
jump distribution in a multi-dimensional setting. 



1. Introduction and Notation 

This article investigates the so called Brownian motion with jump boundary (BMJB) and 
related models, which in recent years gave rise to several interesting results (see e.g. (TUJ, [TT] . 
[Tl] , [B] and [7J). The process itself can be easily described. Consider a Brownian motion (BM) 
with initial value xq in an open domain D C K d which we assume to have a C 2, "-boundary for 
convenience. When hitting the boundary dD of D, the BM gets redistributed in D according to 
the jump distribution v, runs again until it hits the boundary, gets redistributed and repeats this 
behavior forever. Obviously there are a number of possible generalizations, e.g. the BM can be 
replaced by a more general diffusion and the jump distribution might depend on the boundary 
point, which is hit by the diffusion (see [5]). In this work we restrict ourselves to the model of a 
diffusion - corresponding to a 'good' elliptic differential operator L - with a fixed jump distribution 
and partly even to BM with a fixed jump distribution, a setting already leading to interesting and 
highly non-trivial questions concerning e.g. the relation between the spectral gap of the process, 
the spectrum of the Dirichlet Laplacian and the jump distribution. 

The first main result of this paper establishes in a general setting the continuous dependence of 
the spectral gap of the process on the jump distribution, where continuity is meant with respect 
to the weak topology. This answers a question posed by I. Ben-Ari and R. Pinsky in [6] . 

The second main contribution of this work is a direct probabilistic route to the calculation of 
the spectral gap 71 (v) in the case of a one-dimensional Brownian motion in D = (a,b) with an 
arbitrary jump distribution based on the famous probabilistic coupling method. This enables to 
understand the somewhat surprisingly fast rate to equilibrium and its independence on v on a 
'path level'. As was already observed in previous contributions to the theme of the present work 
it turns out that 71 (y) is strictly larger than the lowest eigenvalue Aq '^ of with Dirichlet 

boundary conditions (see e.g. [5] in the general case and [TU], [UJ, [5J and [7J in the case when v 
belongs to the class of Dirac- measures). Even more remarkably the spectral gap 71 [y) is constant 
in v and coincides with the second Dirichlet eigenvalue \^f' b \ This has already been known due 
to a result in [7] in combination with a theorem in [5 . More precisely, Ben-Ari and Pinsky [7] 
had determined the real- valued subset of the spectrum of the generator. Hereupon, Li, Leung and 
Rakesh [5] established that the spectrum is in fact a subset of the real line. Their proof is quite 
involved and uses highly non-trivial results from Fourier analysis. In [5] it is explicitly mentioned 
that the authors can not provide an intuitive explanation for the ^-independence of the spectral 
gap. This remark together with the open question concerning the dependence of the spectral gap 
on the jump distribution - formulated by Ben-Ari and Pinsky in [7] - initiated the present work. 

Despite our purely theoretical study we want to mention that the process considered in this work 
or closely related versions are used in financal applications and also as simple models of processes 
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in neuroscience and operations research. Moreover we point out that some known results together 
with (often slightly different) proofs are included in this work in order to be self-contained and 
readable independent from previous works. 

This work is organized as follows: In section 2, we analyze the L-diffusion with jump boundary 
from a probabilistic point of view. More precisely, we re-prove already-known results concerning 
exponential ergodicity using classical renewal theory, which yields non-tangible upper bounds for 
the convergence rate. A sharp lower bound for this rate is derived in section 3. In section 4 we 
answer the question posed by Iddo Ben-Ari and Ross Pinsky by establishing the continuity of the 
spectral gap 71 (u) in v with respect to the weak topology in a general multi-dimensional setting. 
This continuity property is then used in order to reduce the explicit calculation of 71(f) in the case 
of a one-dimensional BMJB to jump-distributions v with compact support in (a, b). Finally, in 
section 5 we construct an efficient coupling in the case of a one-dimensional BMJB with compactly 
supported jump distribution v. This coupling is used in order to determine 71(f). 

In the remainder of this introduction we fix the notation and recall some basic facts which will 
be needed throughout this work. 

1.1. Notation. We introduce the process in higher dimensions and later specialize to the one 
dimensional situation. Let D C K d be a domain with C 2 ' a -boundary (1 > a > 0) and let 
(f2, (B t )t>o, (^x)xeD) denote a smooth uniformly elliptic diffusion in D which is killed after hitting 
the boundary of D, i.e. (B t )t>o is a diffusion process associated to a generator L of the form 

_^ d d 

L := - ^ a. i j{x)d ij + ^b t {x)d ll 

where the matrix a — (aij)fj =1 is uniformly elliptic with symmetric coefficients 
Moreover, we assume that and 6, are bounded and have bounded derivatives. Of course, 
somewhat weaker assumptions are possible, but we do not aim for such a generality in order 
to avoid technical difficulties. This process induces a compact semigroup of bounded operators 
(Pf)t>o in L 2 (D), which is generated by L, with Dirichlet boundary conditions as operator acting 
in L 2 (D). The spectrum £(— L) thus consists of a sequence (A^)^L converging to infinity. If the 
operator L is symmetric then there exists an associated orthonormal basis (of a suitable L 2 -space) 
( l Pk)T=o °f eigenfunctions. For example L — ^J2i j=i^i( a ij( x )^j) * s symmetric in L 2 (D). We 
denote by g D (-, ■), p D {t 1 •, •) the Green- and the transition-function associated to BM in D killed 
at the boundary dD. If v is a probability measure in D we denote by p D (t, u, y) and g D (v, y) the 
expressions 

/ p D {t,x,y) v(dx) and / g D (x, y) v (dx), 
J d J D 

respectively. Observe that (Pf) t >o acts also in a consistent way on L P (D) for every 1 < p < 00. 

Let W' 1 be a L-diffusion in D with initial distribution p which is killed at dD. Moreover, 

let (W u,i )i>2 denote an independent family of killed L-diffusion in D with initial distribution v , 

which is independent of W* 1 . Set If" = inf{< > | Wf' 1 e dD}, S? = inf{t > | W^ 1 G dD} 

and inductively we define T?^ := Tf'" + S^ +l for i > 1. The process (Af), called a diffusion 

with jump boundary starting from the initial distribution p, is now defined as 

00 

(1.1) xr ■■= i {0 <t<Tr } wf 1 + J2 ^Tt^rryK-Ttr 

i=2 ' 1 

If L is just the Laplacian then we simply write BMJB instead of Brownian motion with jump 
boundary. If the initial distribution is clear from the context we just write (A t !/ ) t >o instead of 

(xr)t>o. 

The process X" induces an operator semi-group on L°°{D) by 

(1.2) VU{x) - E x [/(AT)] = f /(y)P(Ar G dy\X» = x) . 

J a 

Instead of ¥{XP' V G •) we will often write P p (A^ G •) or even P p (A 4 G •)• 
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Moreover, |j • \\tv denotes the total variation norm and for a measure 7 and a measurable 
function / we write occasionally (7,/) instead of J f dj. 

It has been shown by previous authors ([TO], [IT], [T3] and [7]) that this process is uniformly 
geometrically ergodic, i.e. there exists a probability measure /j," such that 

(1.3) -lim ilogsup||P x (xr e-)-^(-)|| T v=7iM>0, 
It is well-known that this is in fact equivalent to 

(1.4) lim log sup \\V!f- [fdn v \\ oo =<n(u), 

and hence 71(f) is referred to as the spectral gap. Since the spectrum £ of an operator and its 
dual coincide, it is reasonable to investigate the spectral properties of V" (t > 0) by studying the 
spectrum of S" on the Banach-space L 1 (Z3), which is often easier to analyze and defined by 

(1.5) L\D) 3g^S»g(-) := / g(x)p v (t,x,-)dx, 



J D 

where p v (t, •, •) denotes the transition density of the process (X^) t >o. The existence of the con- 
tinuous transition density p u (t, •, •) can be deduced from the existence of the continuous transition 
density p D (t, •, •) via decomposing E x [f(X")] according the number of jumps, which occured be- 
fore time t. This is formulated in a precise way in equation (1.2) of 7 . In particular this implies 
that V% maps bounded countinuous function onto bounded continuous functions. 
It can be easily verified that 

(1.6) (S?g,f) = (g,VU) V.g e L\D)J e L°°(D), 
i.e. 

(1.7) sr=v?. 

It is well-known from general spectral theory that this implies 



E(5 t ") = E(P t ") for all * £ R+, 

i.e. up to complex conjugation both operators have the same spectrum at a fixed time t. And 
due to compactness of S% as shown in great generality in [B] one can apply the spectral mapping 
theorem in order to relate the spectrum of the semigroup to the spectrum of its generator (for a 
very readable account on semigroup theory we refer to [8]). 

2. Geometric Ergodicity using Renewal Theory 

In this section we give a short proof of ergodicity and geometric ergodicity of multi-dimensional 
L-diffusion with jump boundary. We want to stress that this result is well-known and not at all 
surprising. We include this result here, since our arguments slightly differ from the ones used in 
previous works and are rather short. The only basic ingredient is renewal theory (see in particular 
and US]). 

We will need the following certainly known Lemma. 

Lemma 1. Assume that we have for every x £ D that T X {X^ E •) — > [i v in the weak topology. 
Then fj, u is the unique invariant measure for X", i.e. 

(2.1) n v (A) = Pp-CX? £ A) for all A e B(D), t e R+. 

Proof. The weak convergence P X (^X^ G •) — > [i v implies that for all bounded and continuous 

feC(D) 

f fdfi = lim E x [/(*£,)] = lim E x [E x [f(X? +s ) \X?\ } 

(2.2) = Km E x [E x ,f(xn] - E„, [/(X t ")] , 

where the last equality follows from the continuity of Ej,[/(X ( !/ )] = J D p"(t,y, z)f(z) dz with 
respect to y e D. One way to establish such this continuity is via equation (|2.3[) . □ 
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Now let us provide a direct route for deriving the form of the invariant measure. 

The above Lemma suggests to calculate E^ [/(X^)] for continuous functions / G C(D). 

E x [f(xn] = Mf( x n,T? > t] +M x [f( x t), n < t\ 

(2.3) =E x [f(X?),T? >t]+J E x [f(X»)\T» = s]F x (T? G da) 



= E x [fW),T? >t]+ [ E v [f(X»_ a )]F g (T? G da) 
Jo 



This is called a delayed renewal equation, which can be further analyzed by considering the 
associated pure renewal equation 

(2.4) E,[f(xn] = Mf(xn,T» >t] + [ K v [f{XZ_ t )]P v (7? G ds), 

Jo 

or for short 

Z(t) = z(t) + Z*F(t), 

with Z(t) = E u [f(X?)], z(t) = E v [f{X%),T? > t] and F(t) = F U (T? G [0,t]). It is well-known 
from general renewal theory (see e.g. [5J) that 

Z(t) = U*z(t), 

where U(t) = YlnLo F* n it) and F*° is defined as the Dirac distribution function degenerated at 
zero. Since z(t) is directly Riemann integrable (for a definition see e.g. [2]), we can apply the Key 
Renewal Theorem (see e.g. [2J, p. 155), which says that 

-1 poo p oo 

(2.5) Z(t) = U*z(t) -> — z(s)ds, m= tdF(t) 



in 



1 f°° 

Km e v [/ro] = y o [/ro, > »] d S 



(2.6) =kW]JJ p D ^v) d ^y) d v- 

f{y) Ejfn dy 

Inserting (|23| into ([23)1 yields 

lim E x [frxnl = / f(y) gD ^ ,y ) dy 

We have just derived the form of the invariant measure in a constructive and probabilistically 
rather direct way. Let us capture the preceding considerations by 

Theorem 1. Let (X^) t >o denote a L-diffusion process in the domain D with an arbitrary jump 
distribution v. Then the process (X^) t >o is ergodic in the sense that for all continuous functions 
f we have 



to obtain 



(2-7) limE4/(Xn] = / Mn v (dy), 

where the invariant distribution fi v is given by 

V"(dy) = -^—g D {v,y)dy, 

with 

(2.8) m»=E v [T?] = f g D (v lV )d y . 

Jd 



ON THE SPECTRAL GAP OF BROWNIAN MOTION WITH JUMP BOUNDARY 



5 



It is natural to ask whether it is possible to obtain more information from renewal theory. At 
this, the speed of convergence to the invariant measure is of main interest. If the convergence 
is geometrically fast with respect to the total-variation norm, the Markov process is called geo- 
metrically ergodic. For general renewal processes, two features of the inter-arrival distribution F 
are necessary and sufficient for geometric ergodicity: Firstly, F must be a spread out distribution 
(for a definition see [2]), since otherwise even convergence in total variation fails to be true. Sec- 
ondly, F must have exponentially decreasing tails (see e.g. [5]). In our setting the inter-arrival 
distribution F V (T^ € •) share both properties, since F V (T" € •) even admits a density. 

Hence [2, Theorem 2.10] allows to refine (|2.5|) to the assertion 

1 f°° 

(2.9) Z(t) = U*z(t) = — z{s)ds + 0(e~ et ) 



for some sufficiently small e. In this situation, the function / in the definition of z(t) is only 
required to be bounded and measurable. Therefore, equation (|2.9[) implies geometric ergodicity of 
the Markov process X t . Actually, the tails of F X \T" 6 •) are uniform-exponentially decreasing for 
x G (a, b), and hence uniform ergodicity can be deduced. 

Despite these affirmative results, it is impossible to obtain good or even sharp convergence rates 
from general renewal theory directly, since the above e has been deduced from a general coupling 
argument which does not fit to our particular situation. 

Similar calculations to the ones presented above have already been used in previous works of 
Grigorescu and Kang pTO] as well as Ben-Ari and Pinsky [7J in order to show ergodicity of the 
process. Grigorescu and Kang [TU] consider the Laplace transform of an expression similar to (12.31) 
and calculate the associated poles. As one expects from general operator theoretic principles the 
density of the equilibrium measure can be characterized as residuum of the pole at 0. Ben-Ari 
and Pinsky use a different but still quite analytic approach. A probabilistic approach using the 
Doeblin condition in a general setting is presented by Grigorescu and Kang in |12j . 

3. AN UPPER BOUND FOR THE RATE OF CONVERGENCE 

The next natural question consists of calculating the exponential convergence rate. Let us look 
at the one-dimensional case. In order to find the exponential rate of convergence the analyst 
usually works on the generator level and considers the following eigenvalue problem 

1 . . f b 

(3.1) u = Xu with 'boundary' condition u(a) — / u(y) v{dy) — u(b). 

2 J a 

The exponential rate of convergence 71(f) coincides with min{3?A 7^ | 3u\ ^ satisfying (I3.1[) }. 
Observe that (13. 1 [) includes the non-local boundary condition, which represents the jump mech- 
anism of the process. Instead of investigating the stationary problem (|3.1I) . we work directly on 
the path level. In this section we first prove a simple upper bound on the convergence rate in the 
case of the one-dimensional BMJB in an interval D = (a,b). For this reason we introduce the 
reflection map R defined by 

R : (a, b) — ¥ (a, b), R(x) := a + b — x. 

Furthermore recall from section 1.1 that (a!l )£1 denotes the sequence of eigenvalues of the 
operator in (a,b) with Dirichlet boundary conditions at both endpoints. 

Proposition 1. Let d = 1 and let (X t ) t >o denote the BMJB in (a, b) with jump distribution v. 
Then we have 

7iH<A^ } . 

Proof. Due to the trivial inequality 

(3.2) \\F x (X t G-)-¥ R(x) (X t e-)\\ TV <2 sup ||P x (JT t G-)-A*|lrv 

it is enough to estimate the left side of (|3.2p from below. Observe first that by symmetry we have 

(3.3) F X (T? e •) =W R{x) (T? e •)• 
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Using (|3.3I) in the last step we obtain 

F x (X t G A) - F R(x) {X t G A) = F x (X t G A; T{ > t) - F R{x) (X t G A; T" > t) 

+ f P„(X t _ r G A)P x (T 1 i/ Edr)- [ P„(X t _ s G A)P R(x) (7? G da) 

Jo Jo 
= P x (B t G A; T (Q)b) > <) - P fi(:c) (Si G A; T (o , 6) > i), 

where (B t )t>o is a one-dimensional Brownian motion and T( a ^) = fnf{i > | B t (a, b)}. From 

the spectral decomposition of the killed BM in (a, b), <p^ (x) = an d ^2i+i( a; ) = 

-V^+l(i2(a;)) for alH G N U {0} we deduce 

F x (X t E A) — F R(x) (X t e A) = F x (B t G A; T (o>6) > t) - P fl(x) (S t G A; T (aM) > t) 

z=0 JA i=0 

OO „ 

i=o Ja 

This together with (13.21) clearly proves the assertion of the proposition. □ 

One should keep in mind the important role of the symmetry (x,R(x)) used heavily in Propo- 
sition (|T|). This symmetry will occur also later in the proof of the matching upper bound. 

Thus the straightforward arguments used in this section demonstrate in a rather direct way 
that in the case of a one-dimensional BMJB with jump distribution v 

(3.4) 7iH<A^ } . 

Before we start with a more precise investigation of the spectral gap for general jump distributions 
we ask whether the upper bound in (|3.4[) is sharp at least in a very special situation. We have 
seen that the L-diffusion with jump distribution v is ergodic and the unique invariant distribution 
\i v has the form 



m Jd 



'D 

where m" = j D J D g D (x, y) v(dx) dy. It seems reasonable that the speed of convergence to equilib- 
rium should be fast whenever jump- and invariant distribution coincide. This immediately suggest 
the question for the existence of a jump measure such that v = [i v '. For simplicity let L now be a 
symmetric operator, for which the theory of quasi-limiting distributions is more transparent. Still 
the result also holds up to obvious modifications in the non-symmetric case, but the corresponding 
results on the quasilimiting behavior are scattered widely in the literature. It turns out that there 
is a suitable jump-measure, namely the quasistationary distribution for the killed L-diffusion in D, 
which is given by po '■= f D (y) dy (normalized to give total mass 1). The notation quasi-stationary 
refers to the fact that for measurable A C D 



[W t e A | T D > t) = f <ffi{y) dy = p D (A) 

J A 



A 

where (Wt)t>o denotes an L-diffusion and To = inf{t > | Wt G 3D}. This fact has already 
been noticed by Ben-Ari and Pinsky in Theorem 1 of [5], where it has also been shown that the 
quasi-stationary distribution defines the unique fixed point of the equation fi v = v. 

Hence if (|3.4p would be tight, it should be at least tight for the BMJB in D with jump distri- 
bution v = pi)- That this is actually the case is the content of the following 

Proposition 2 ([B], Theorem 1). Let L be a symmetric operator in D and let the jump distribution 
v be given by the quasi- stationary distribution pn(dy) — tp D (y)dy. Then we have 

7iM=Af. 
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Proof. The proof differs from the one given in [6] and seems to be more probabilistic to us. Let 
(X t )t>o denote the L-diffusion in D with jump distribution v(dy) = ipo(y)dy. For measurable A 
we decompose F x (X t 6 ^4) according to whether there has been a jump before time t: 

V x (X t G A) = F x (X t G A. T" <t)+ F x (X t e A, T» > t) 

= / F v (X t _ s G A)F X (T? G ds) + P x (X t G A | Tl > f)P a (2? > t) 
Jo 

= v{A) - v{A)F x (T" >t)+ P x (X t G A | T" > t)F x {T? > t) 
= v(A) + (V x (X t G A | Tl > t) - v{Aj)F x {T? > t). 

Thus observing that - lim^^ \ log \\F X {X% G • | T% > t) - v{-)\\tv = Af - Xq (see e.g. [16] and 
[5]) we arrive at the assertion of the Lemma. □ 

4. Continuity of 71 (v) 

Rather simple symmetry-arguments have been used to prove that 71 (v) < X\. For the exact 
evaluation of 71 {v) a deeper analysis is necessary. One natural approach consists of reducing the 
problem to 'simple' jump distributions by certain continuity considerations. The fact that the 
equilibrium measure \x v depends continuously on v with respect to the weak topology raises the 
question, whether the same is true for the spectral gap A(V). The answer is affirmative and we 
want to stress that the results of this chapter are valid in the multi-dimensional situation. 

Instead of dealing with V" directly, we investigate the behavior of the semigroup S" defined in 
(|1.5|) , where we know that S^* = V\ ■ The main reason for this procedure is that in contrast to 
V" , S" turns out to be strongly continuous (see [5] and [7]) so that the general theory of strongly 
continuous semigroups (see e.g. [8]) can be applied. For example, this makes it possible to verify 
compactness for and to relate its spectrum to the spectrum of its generator by application of 
known results of semigroup theory (see e.g. [7]). The following two basic results from operator- 
theoretic perturbation theory, which are due to Kato, will be applied later on to the semigroup 

Theorem 2 f[13j. p. 178). Let T be a closed operator in the Banach space X and let o~{T) be 
separated into two parts S'(7~) and £"(T) by a closed curve T, i.e. a rectifiable, simple closed 
curve r can be drawn so as to enclose an open set containing S'(7~) in its interior and £"(T) in 
its exterior. Then one has a decomposition of T according to the decomposition of X = M'{T) © 
M"(T) of the space in such a way that the spectra of the parts Tm' md Tm" coincide with £'(7") 
and £"(7~), respectively and 7m' £ B{X). 

In the following <5(-, •) denotes the metric which defines the generalized convergence of closed 
operators (compare [T5] p. 202). 

Theorem 3 f[13j. p. 212). Let T be a closed operator in the Banach space X and let o~(T) be 
separated into two parts S'(T) and X"(T) by a closed curve V, i.e. a rectifiable, simple closed 
curve r can be drawn so as to enclose an open set containing S'(7~) in its interior and S"(T) in its 
exterior. Let X = M'(T) ®M"(T) be the associated composition of X. Then there exists a S > 0, 
depending on T and T with the following properties. Any closed operator T on X with 5(T, T) < S 
has spectrum tr(T) likewise separated by T into two parts Y/(J~) and £"(7~). In the associated 
decomposition X = M'{T) ©M"(T), M'(J~) and M"(T) are respectively isomorphic with M'{T) 
and M"(T). In particular dim Af '(f) = dim M'(T) and dim M"{f) = dim M"(T) and both 
£(T) and S"(7~) are non-empty if this is true for T ■ The decomposition X = M'(7~) © M"(T) 
is continuous in T in the sense that the projection P[T] of X onto M'(T) along M"(T) tends to 
P[T] in norm as 5(7", T) — > 0. 

In order to apply Theorem [3] observe that for a bounded operator T and a sequence of closed 
operators (7~ n ) n we have by Theorem 2.23 in |13j that lim n ^ x S(T n ,T) = if and only if the 
operators T n are bounded for n large enough and lim^oo \\T n — T\\ = 0. Now, Theorems [5] and [3] 
immediately imply the following conclusion. 
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Corollary 1. Assume that (T n ) n <£N is a sequence of compact operators in a Banach space X 
which converges with respect to the operator norm to a (necessarily compact) operator T and 
let r be a curve such that T divides the spectrum £(7~) of T into a two parts, one finite part 
Y,'{T) lying inside the curve with no limit point and one part £"(T) lying outside this curve. 
Then for n large enough the spectra S(7^) of the operators J~ n are also separated by V and since 
dim AI'{Tn) = dim M'(T) < oo and since T n P[Tn\ converges in operator norm to TP[T] as n — > oo 
we also have S'(T n ) — > £'(T) as n oo with respect to the Hausdorff metric and in the following 
sense: 

Ve > 03N e NVA G E'(T)Vn > iV 3A?, . . . , A£ £E'(T»): max |A - Af| < s. 
Furthermore, If m (m™ ) denotes the algebraic multiplicity of the eigenvalue A (Xf) ofT (T n ), then 

Ef=i m ? = m - 

Now we turn to the application of these abstract results to the situation at hand. We want 
to stress that several of the following arguments are taken from [TJ. We give rather detailed 
arguments since some care concerning the dependence of various estimates on the jump distribution 
is necessary. 

We will need the following Lemmata. These are probably well-known, but we have not been 
able to locate them in the literature in the form we use them. 

Lemma 2. Let D C K d be a bounded domain and let as above (A^ ) ) n gN an d ('Pn )ngN denote the 
sequence of eigenvalues (counted according to their multiplicities) and the associated eigenf unctions 
of the symmetric operator L 



1 

L ■= ^ Yl d i( a ij(-) d j) 
in D with Dirichlet boundary conditions, respectively. Then for any e > we have 



(4-1) ^e-^||^|| 2 00 <oo. 

k=0 

Proof. By Jensen's inequality we have 

(^(x)) 2 = e 2 ^^ x [^{B t ),r D >t}) 2 < e 2 ^E x [^(B t ) 2 ,r D > t}. 
Since swp x y&D p D {t, x, y) < t~ d l 2 (see e.g. pQ), we conclude that 

(A o\ \\,J3\\1 „2A fe i + — d/2 ||,„,D||2 „2\ k t ,-d/2 

Choosing t = | in (|4.2|) . the finiteness of the sum in (|4.ip follows from Weyl's theorem, which states 
that for appropriate constants c\,C2 > one has c\k 2 l d < \® < C2k 2 / d (see e.g. [3], Theorem 
6.3.1). □ 

A simple application of Lemma [5] and the dominated convergence theorem gives 

Corollary 2. Let D C M. d be a bounded domain with smooth boundary and L a symmetric diffusion 
operator with Dirichlet boundary conditions. Let Trj = mf{t > | B t G dD} the first hitting time 
of the boundary of D the L-diffusion process (B t )t>o starting in x G D then the distribution of To 
is absolutely continuous with a density which is jointly continuous in {x,t) G D X (0,oo). 

Let Lq = |Ei j=odi( a ijdj) ^> e t ne generator of a uniformly elliptic reversible diffusion in M. d 
with smooth and bounded coefficients written in divergence form or 

d -, d 



i,j—l i—0 i—0 

in non-divergence form. We have seen that the first hitting time distribution of the set dD has 
a density which is continuous with respect to the starting point of the Ln-diffusion. We will now 
use this result in order to extend this to diffusions of the form 

L:=Lo + 7'V, 
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with 7 e (C c {R d )) d . 

Lemma 3. Let V x denote the law of a L-diffusion and let as above Td = inf{t > | Bt G dD} 
denote the first hitting time of dD. Then T x (Td G •) is absolutely continuous with a density h{-, •) 
that is jointly continuous in (x, t) G D x (0, oo). 

Proof. Let us denote by V x the path space measures associated to the diffusion generated by L and 
let P x denote the path space measures corresponding to the Lo-diffusion. Then by the Girsanov 
theorem (see e.g. section 1.9 in [IT]) one has 



dP° x 
where 



1 '* 



7V t :=exp( / a- 1 1 {X s )dX s -- / (a^TpQ, 7 PQV ds 



where X t = X t — J * P(X S ) ds. Therefore we have 

Px(?fl>i) = E°[iV t ;T D > t], 

since we already know that the exit time distribution P°(Td 6 •) of the reversible diffusion 
corresponding to L is absolute continuous with respect to the Lebesgue measure the same holds 
true for P x (Tq G •). Thus there is a function h(x,t) (x G D,t > 0), which is jointly measurable 
such that for x G D and measurable A C (0, oo) 



: (Td e A) = / h{x,s)ds. 

J A 



It remains to show that h(x,t) has a continuous modification in (x,t) E D x (0,oo). As a first 
step let us show that h(x,t) is a.s. uniformly bounded in x for all t > c > 0, i.e. h(x,t) < M t < 
M c < oo Vt G [c, oo]. We will use Lebesgue's version of the fundamental theorem of calculus to 
conclude that for almost every (x, t) G D x (e, oo) 

(4.3) lim 7T T5T m / P *( T £» G (t - u,t + u)) dx = h(x,t) 

r,u^o2u\B(x,r)\ J B(x r) 

as well as 



1 1 

^]B(x~r) 



P x (T D e(t-u,t + u)) dx 

B(x,r) 

1 1 



2u \B{x,r)\ 



D 

V 

B{x,r) J D 
11/ IV, 



p (t — e,x,y)P y (TD G (e — u,e + u)) dy dx 



(4.4) =TTTW \\ p D (t-e,x,y)P y (T D e(e-u,e + u))dydx 

ZU \&{X, rj| JB(x,r) 



D 



1 1 



2u \B(x,r) 

D 



e+u r r- 

/ / p D {t - e,x,y)h(y,s)dydxds 

e-u J B{x,r) J D 



as r, u — y 0+. Using equations (|4.3I) and (|4.4j) together with the fact that for very positive t > 
the operaror Pf maps an integrable function to a bounded one we deduce that for some constant 
C > 



\h(x, t) — h(z, s)\ = 



{p D {t - u,x,y) ~p D {s - u,z,y))h(y,u)dy 



(4.5) < C \p D (t-u,x,y)-p D (s-u,z,y)\dy 



/' 

£> 

1^ 

£> 
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Due to continuity of p D (-, ■, y), we see that for (z, s) — > (x, t) G D x (0, oo) the integrand in (14.51) 
converges pointwise to zero, and hence we obtain by Lebesgue's theorem that 

lim \h(x, t) — h(z, s)\ = 0, 
;z,s) >(x,t) 

which yields the claim. 

□ 

Lemma 4. Let (v n )n<£N a sequence of probability measures in D converging weakly to the probability 
measure v. Then for every e > 

lim sup \p D (s, v n ,y) -p D (s,v,y)\ = 

n ->°°s>E,yeD 



Proof. First observe that due to limn^oo v n = v weakly we have for every s > and z G D 

lim p D (s, v n , z) = p D (s, v, z) 

n— Yoo 

as well as 

lim / p D (s,u n , z) dz — I p D '(s,v, z) dz. 
™^°°Jd Jd 
Thus as a consequence of Scheffe's lemma we conclude that 

(4.6) lim / \p D (s,v n ,z) — p D (s,y,z)\dz = 0. 

n^ocJ D l 

Observe now that due to the Chapman-Kolmogorov equations 



p D (e + s, v ni z) -p u {e + s, u, z) \ = 



(p D {e/2, v n , y) - P D (e/2, v, y))p D {e/2 + s, y, z) dy 

D 



<[ sup p D (e + s,y,z)] / \p D (s/2, v n ,y) ~p D {e/2, v, y)\dy. 

s>0;y,zeD Jd 

Since according to equation (|4.6p the right hand side converges to the assertion of the Lemma 
is shown. □ 

Lemma 5. Assume that i>i — > v. Then there exists for each t > an a — a(t) < 1 such that for 
all v G {u, Vi,i>2, ■ ■ ■} an d x G D 

P^C^+i < t) < a 1 Vi G N. 
Proof. For v G {v, vi 1 i>2, ■ ■ ■ } define 

c-v = Cu(t) = 1 - / / p D (t,x,y)dyv(dx). 
Jd Jd 

As a direct consequence of the Markov property we get for i> G {y, v\, 1/2, ■ ■ ■ } 

P*(lti < t) < (P(3 < a < < = W?' 1 G 
= 4 

^From limn^oo v n = v it follows that for all i > we have 

Cj/ n — > c„ for n — > 00, 

which implies that for all e > there exists n(e) with the property that for n > n(e) 
(4.8) c Vn < c v + e. 

For e small enough (e.g. e < 1 ~ c " ) the right hand side in (I4.8[) is strictly smaller than one. Since 

c Vn < 1 for all ra € {1,2, ... , n(e)}. 

Let 

Obviously, a — a(t) < 1 and since c v < a, the result follows from (14.71) . □ 
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D JD JO 



Before we state and prove the main result of this section, let us introduce important notations, 
which are necessary for illustrating the structure of the proof. Let 

(4.9) h"(t) = |P P (T» < t) 
and for n > 2 

(4.10) K»{t) = |p(2r < t) = (h" * (ft")*'- 1 )^). 

In the case when p = 5 X is the Dirac measure in x, we will write h% v {t) instead of h s ^^(t). Observe 
that h p (t) is nothing else than the density of the first exit time of the L-diffusion with the initial 
distribution p from the domain D. Moreover, set for n > 1 

s o.t9( z ) '■= I p D {t,x,z)g(x)dx and S^ t g{z) := [ [ p D (t - s,y, z)g(y)h% u (s) ds v{dy)dx 
and similarly 

VZJ(x) := f p D (t,x,y)f(y)dy and T v nt f{x) := f f p D {t - s,y, z)f(z)h x n > v (s) dsv(dy)dz. 

JD JD JO 

With these notations let us recall the definition of P" and write it in another way. For all / G L° 
we have 

oo 

vtf(x) = e x [f(xn] = e x [f(xn,T? >t\+Y,®* [f(xn, n<t< r* +1 ] 

71=1 

oo „t 

= E x [f(X?),T? >t]+J2 Mf( X t),K +1 > t\T» = s]P x (TZ G ds) 

n=l Jo 

oo -t 

(4.11) =E x [f(xnX >t] +Y, / Mf( x ts);S^ +1 >t~ s ]h^( s )ds 

n=l ^ 

= / p D (i,x,y)/( 2/ )dy + V /"*/"/ p D {t- s,y,z)f(z)v{dy)dzhY{s)ds 
Jd n=1 Jo in J_d 

oo 

and hence 

oo 

(4.12) P"=X>m- 

i=0 

A dual version of this calculation immediately yields 

oo 

(4-13) S» = Y. S tf 

Moreover, it is obvious that 

(4-14) ($«>,/) = (<?,7^/) V ff G L 1 , / G i G NU {0}, 

and hence 

«S£ = VzeNu{0}. 



We are now able to prove our first main result. During the proof we use || • \\ PtP to denote 
the operator norm of a bounded linear operator mapping L p to L p . 

Theorem 4. The spectral gap 71(f) of the generator of the semigroup (<S f ")t>o acting in L l (D) 
depends continuously on the jump distribution v with respect to the weak topology. 

Several ideas of the proof are taken from [7|, where they are used for a different purpose. 
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Proof. We will use Kato's results cited in the beginning of this section. In order to apply these 
results we have to show that for admissable distributions v, v i, v%, . . . 

lim v n = v with respect to the weak topology => lim \\S" n — <S^||i,i = 0. for all t £ M+ 

n— too n— >oc 

We have 



N 



N 



N 



- Hi,! < \\St - E^lli.i + IIE( 5 M ~ 5 m)Hm + W S " - E 5 M Hi,! • 



i=0 



i=0 



i=0 



I(N,i>) 

By duality and Lemma [5] we obtain 



II(N,v,u„) 



III(N,u n ) 



(4.15) 



W") = II E 5 mIIm = II E ntL.oo 

i=JV+l t=iV+l 

sup E a (/(X t ");T£ +1 <i) 

/ei°°:||/IU=i 

< supP x (T£ +1 <t) ^aft)*. 



Observe that by Lemma O inequality (|4. 1 5[) holds for all v £ {i>, v\, v<z, . . . , i> n , . . .} with v n — > v 
(in the weak topology). Hence, for all e > we can find no = no(e) such that for N > no we have 



(4.16) 



I{N,v) < e, III[N,v n ) < e. 



Since II(N, v, v n ) < ArrnaXj e {i ) 2 1 ... )J io}{|| < Sfi — <STt Hi i}> ^ remains to show that for all e > there 
exists m such that 

max AWKt ~ S id II 1,1 > < T7 Vn > m. 



ie{l,2,...,n } 



AT 



Since we maximize over finitely many elements, it is easy to see that it suffices to show that every 
e > and every i$ G {1, 2, . . .} there exists ni such that 



(4.17) 

We have by duality 



\\Sl t -S£ t \\ ltl <-Vn> ni . 



\K,t-K:th,i = \\K,t-K: t \Uoo 



sup sup 

ajel3/Gi c »:||/|| o=l 



p D (t - s, ^, z)/(z) dyh^(s) ds 



'0 •/£> 



- s,v ni y)f{y)dzh x ^ Vn (s) ds 



io Jd 
ft 



< sup sup 

xED f£L°°:\\f W^^l JO JD 



sup sup 

xeB/eL°°:||/||oo=i 



p D (t - s, i/, z)f(y) dzti*' v (s) ds 



p D (t - s, v n , z)f(z) dzh^(s) ds 



o Jd 



p (t - s,v n ,z)til' v dsf(z)dz 



o Jd 



p D (t - s, u n , z)ti?' Vn (s) dsf(z) dz 



d Jo 



sup sup 

se£>/ei°°:||/|U=i 



(p D (t - s,u,z) -p u (t - s,v n ,z))f{z)dz I h x ' u {s)ds 



\JD 

/ 



+ sup sup 

xeD/eL«":||/||oo=i 

/„(/)+ !/„(/). 



J D 



j? D (i - s,^„,z)/(z)dz(^ o ,,y (s) - ^'""(s)) cis 
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The first term can be estimated by 



In(f) < SUP 

xeD 



= sup 

x£D 



\JD 



\p (t-e,u,z) -p (t- s,v n ,z)\dz\ h^ u (s)ds 



a(n,t — s) 



< sup 

x£D 



a(n, t — s)h^' v (s) ds 



a(n, t — s)/i io ' (s) ds\ + sup 



a(n, t — s)h^' v (s) ds 



where e £ (0, t) is arbitrary. Using a(n, t — s) < 2 for every n € N and J*h% v {s)ds < 1 for every 
x 6 D we conclude that 



(4.18) 



At(/) < SU P a ( n ) i — s) + 2 sup / l\^(s) ds. 

S<t—£ xGDJt — E 



Due to Lemma H] the first term in (|4.18[) converges to for every e > 0. Whereas the second term 
in (|4.18p can be made arbitrary small by chosing e small enough (see Lemma 3 in [7]). Thus 

limsup /„(/) = 0. 

n— >oo 

In order to estimate //„(/), let 

(4.19) ~h% = h x * {h v * ■ ■ • * h v ) * (h u " *• •.*/»"») . 

„ ' y v •* 

j terms ia-j— l terms 

^From commutativity of the convolution it follows 

■X 



ft? -/^_i iio _j(s) = ft x * * ■ ■ ■ * /iX fe" - * ft"" * • • • * h""] {s) 



(4.20) 

Therefore we have 



j-i terms i -j-i terms 

(h v - h u ") * h x * [h" * ■ • • * h" * h Un * ■ ■ ■ * h u "]{s) 

(h»-h»")*~hU,io-j-l- 



SU P / \ h io U ( s )- h io Vn i s )\ ds 



Iln(f) < 

< sup^ / l/if+i^^-a-^io-i-ilWda 

i6D „•_ n JO 



x£D Jo 

if) -2 





sup 




xED 


< 


(io ~ 


Fubini 




< 


{io - 


< 


{io - 



3=0 
io— 2 ~t 

^ 

3=0 J ° 



lX j,i a -j-2 
I 



(s) ds 



eD3e{o,i,...,i -2} J 
t 



\(h v -h v »)\*hl i0 _ j _ i (s)ds 



x£D Jo 



(4.21) 



— » for n — > oo, 



where jo denotes the index where the maximum is attained. For the last step we use that according 
to Lemma [3] the density 

x (T D e ds) 



h(x, s) 



ds 
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of the first exit time To from the domain D is continuous in x, and hence by weak convergence 
the integrand of the inner integral converges pointwise to zero. Another application of Scheffe's 
lemma shows that the last integral in (|4.2ip converges to zero. From (|4.18p and (|4.21[) it follows 
that there exists m such that (14.171) holds. □ 

Remark 5. 

(1) In [5J Ben-Ari and Pinsky formulated the question, whether the spectral gap 71(f) depends 
continuously on v. Theorem^ answers this question affirmatively. But observe that The- 
orem [7] actually shows that arbitrary finite subsets counted according to multiplicities of 
the spectrum depend continuously on v. 

(2) In the recent preprint [15] we also give answers to Question 1 and Question 2, which are 
posed on page 130 [5]. 

The continuity of the spectral gap might be of some importance e.g. in optimization problems 
concerning the speed of convergence for the multi-dimensional BMJB. 

Corollary 3. Let D C M. d satisfy our standard assumptions and let 71(f) denote the spectral 
gap of the BMJB in D with jump distribution v . Then for every precompact (with respect to the 
weak topology) subset K of probability measures on D satisfying sup nelc 71(f) < 00 there exists a 
jump-distribution vq such that 

71(2/0) = sup 71(2/). 
ueK 

Proof. Take a sequence {v n ) n C K, such that limn^oo 7i(2/ n ) = sup„ K 71(f) < 00. Then by 
precompactness of JC we can extract a weakly convergent subsequence (v nk )kefi with limit fo and 
due to the continuity of the spectral gap 

71(^0) = I™ li(Vn k ) = sup 71(f). 

□ 

For the remaining part of this work we note the following simple 

Corollary 4. Let D = (a, 6) and let v be an admissible jump distribution, i.e. f({a, b}) = and 
set f„(-) = f(- I {x ; dist(x,dD) > l/n}). Then we have 

lim 71 (f„) = 71(2/). 

n— >-oo 

Thus once it is shown, that 7i(-) is constant on compactly supported jump distributions (or 
even on jump distributions supported on finite sets), we can conclude that 7i(-) is constant on the 
set of all admissible jump distributions. 

5. Rate of Convergence: Probabilistic Approach 

In this section we consider BMJB on the interval (a, b) and recover a recent result of [5], 
which was previously shown via elegant Fourier-analytic arguments. As already mentioned in [S] 
these Fourier-analytic arguments do not offer any probabilistic explanation. We use the coupling 
approach, which provides a probabilistically more satisfactory explanation for the obtained con- 
vergence rates. As is well-known the coupling method usually involves the construction of two 
suitably dependent processes as a crucial step. In our situation the main difficulty is present due to 
the fact that in contrast to one-dimensional diffusions without jumps two independent BMJB can 
'pass' each other without hitting each other. The lower bound on the convergence rate is derived 
- as usual in coupling approaches - via an investigation of the tail behavior of the coupling time. 
The following result constitutes the second main theorem of this work. 

Theorem 6. Assume that d = 1. Then we have 

(1) 7iM = Af 

(2) There exists an efficient coupling. 
The proof is split into smaller pieces. 
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5.1. Auxiliary results. We will need the following elementary auxiliary results: 

Proposition 3. Let I — (a, b) be an open interval with center c — For y £ I let r y — inf{i : 
y + B t £ dl} the first time of leaving the interval. Then we have for all y £ I and t £ K+ 

P(r y > t) < P(r c > t). 

Proof. The proof is based on a simple coupling argument: Without loss of generality we may 
assume that b > y > c. We define the coupling as follows: For t < r cy = inf{i : X t — Y t } let 
X t — y — B t and Y t — c + B t . Now let us distinguish the following two cases: First case: X and 
Y meet in at time r cy < t c . In this case let Y = X = + B t — B T i for t > r cy and hence 
both processes leave the interval at the same time. 

Second case: X and Y do not meet each other before t c . In this case, we have by definition of the 
processes that r y < t c . Hence we obtain 

P(r c >i) = V(T c >t,T c >T Cy )+V(T c >t,T c <T C y) 
= P(Ty >t,T c > T Cy ) + P(T C > t, T C < T Cy ) 
> P(r y >t,T c > T Cy ) + P( Ty >t,T c < T Cy ) 

= P(r„>i). 

□ 

Now we come back to the symmetry argument, which already played an essential idea in Propo- 
sition [TJ In the following Proposition we put this symmetry argument - which is the key idea in 
our approach - in the coupling context. 

Proposition 4. Suppose that a < x < c = < y < b and y — R{x) = a + b — x. Then there 
is a coupling of X u,x and X u,v such that the coupling time is equal to the exit time of a standard 
BM (=starting at 0) from the interval I = I(x, y) = {—^-^-, ^ b ~°^~^- x ~ y ^ ) . 

Proof. We use the same notation as in the previous proof and set r = t x A r xy . For t < r define 
X v - X = x + B t and X v,v = y — B t . Since y — R(x), we have to distinguish 



• First Case: 

X V T - X = x + B T = y-B T = Xp y = c. For t > r we define X v t x = X"' v = c + B t _ T . 

• Second Case: 

x + B T = a and y — B T = b. For t > t we define X"' x = X^' v — v + B t - T 
We see that as long as the path of B t is contained in the interval / = (— ' b ~ a )~( x ~ v ' ), the 
processes X v,x and X v > y do not merge, but once the Brownian motion (B t )t>o exits dl they 
immediately colesce. This yields the claim. □ 

Remark 7. A crucial observation is that \ J(x, R(x))\ = is independent of x and that = 

Now we will show that this idea can be extended to processes with arbitrary initial values x, y 
whenever the jump measure is compactly supported in (a, 6). 

5.2. Construction of the coupling. Now let us built up the coupling. First of all note that by 
a very simple argument using the triangle inequality and a symmetry argument we can assume 
without loss of generality that 

cl ~\~ b x ~\~ y 

(5.1) < y — x < dist(supp(v), {a, b}) and — - — < — - — . 

Let us introduce two copies of BMJB, called X and Y in the following way: 

a) Let xi = x, y% = y and X t = x\ + B t , Y t = yi + B t . Now stop stage a) at time 
n = r sym A n, where r sym = inf{i : R(xi + B t ) = yi + B t } and n = inf{i :y 1 + B t = b}, 
i.e. we stop when either the copies are in a symmetric position (t\ = r sym ) or when Y hits 
b (ti = Tf,). Note that due to the assumption (|5.1|) one of these two cases has to occur. 



16 



MARTIN KOLB AND ACHIM WUBKER 



First case: If t\ = T sym then at time t\ we are in the situation to apply the coupling 
presented in Proposition 01 

Second case: If t\ = t\, then we have x 2 := X Tl = X Tl - = x\ + (b — y) = b— (y — x) and 
Y Tl — 2/2 = Jii where J\ denotes the first jump with distribution v. The construction is 
continued with the next stage, 
b) Let X = x 2 + (B t - B T1 )(= x\ + B t ) and 

(5.2) Y = y 2 — (B t — B T1 ) = f (-ljWodB.. 



n 



To see that the stochastic integral in the definition of Y is well-defined, let Tt = cr(B s , s < 
t) and Qt = cr(Tt,cr(Ji)) the cr-field generated by er(Ji) and Tt- Now observe that 
(— 1) 1 {ti<*} is Q t -measurable and that by independence of er(Ji) and Tt B t is still a Brow- 
nian motion with respect to Qt, the integral in (I5.2j) is still well-defined. Now recall that 
X2 — 2/2 = x + (b — y) — Ji = b — (Ji + (y — x)) > 0. Stop stage at time t 2 = either copies 
meet or when the distance from Y to X is b — J\ . 
First case: The coupling occurs and we are done. 

Second case: The initial distance from Y to X is (b— Ji) — (y — x). Therefore, Y moved 
^ to the left, i.e. 2/3 = 2/2- ^ = Ji - ^ > a (due to (jOl ). and X moves ^ to 
the right, i.e. X stops at x^ — x 2 + iL ^- — b — < b. Now move to stage 3. 
c) Let X — x 3 + (B t — B T2 ), Y = y 3 + (B t — B T2 ). Recall that £3 > y 3 and x% — 2/3 =b — J\. 
Stop stage c) at time T3 = either X and Y are symmetric or X hits b. 
First case In the symmetric case we continue the coupling as in Proposition 2] 
Second case: Here we have X T3 - = b and Y T3 = b — (6 — J\) = J\. At time T3 we let X 
jump to J\, i.e. we use the same J\ as has earlier been used for Y. Hence the coupling 
occurs and we are done. 

First let us observe that the above construction is in fact a coupling for the processes, i.e. both 
processes have the same marginal distributions. Denote Thu = vaf{t : X t hits the boundary dl}. 
Then, for t < Thit we have by construction of X t 

x + B u 

since all changes of sign of the BM appear by definition in the F-process. This immediately 
implies that X, up to its first arrival to the boundary, is independent of J\ and behaves as a BM. 
It remains to show that Y also behaves as a BM in the interior of the domain. Now let 77^ (r^ it ) 
be the first (second) time when Y hits d(I) and observe that by construction for t < rh t , Y can 
be written as 



and similarly for t% u > t > r\. 



hit 



Y t =y+ [ f(s)dB s 
Jo 

Y t = Ji + [ f{s)dB s 



where f(s) is C/ s -measurable and only takes values in { — I, I}, because / keeps track of the changes 
of sign of the BM as defined in the coupling construction. But up to time t these changes only 
depend on the Jj-history of the BM and on a{J\). But from here it is easy to see, e.g. using 
Levy's characterization of the BM that Y also behaves as a BM in /. 

5.3. Analysis of the construction. The crucial observation in the analysis of the coupling 
construction is that in each stage the stopping rule is determined by an exit time of a standard 
BM on an interval of length less or equal to . These exit times (which are not independent) 
can be dominated by a sum of independent exit times of BM's of an interval of length as we 
will see in the sequel. 

a) First case: In this case we stop the process when the two copies are symmetric about 
the origin. For this to occur the BM has to arrive to ( a + fc )~fa~ x ) _ x (terminal — initial 
initial location of X)= (a+b) ~ fa+z) < 0. 
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Second case: Here we stop when the coordinate Y hits 6, i.e. in the case when the BM 
B has hit b - y > 0. 

(y+x)-(a+b) b-a 



Thus, stage a) ends when the BM exists an interval of length b — y + {y+x> {a+ ' < 
b) First case: We stop when X and Y arrive at the point — X2 ~ V2 — (y- x )-( b - J i) 



2 



Second case: We stop when BM arrives to 2!— H > Q 



2 2 

__ -. II 



Thus, stage b) ends when BM exists an interval of length b 2 Jl < 

c) First case: We stop when X and Y are in symmetric locations. Initially, X is in 
X3 = b — and X is kept b — J above Y throughout the stage. Therefore we stop 
when the BM arrives to (a+b) + (b ' J) - x 3 = a+ ^~ J < 0. 
Second case: We stop when BM hits b — x% — 



2 



Thus, stage c) ends when BM exists an interval of length < . 
Now we are able to finish the proof of Theorem [5] by estimating the tails of the coupling time. 

Proof. (Theorem [5]) Now let £i, £21 • • • i £5 De independent exit times of standard Brownian motion 
of the interval dl = (— ^ip, 9 ^). Then, by Proposition [3] and Proposition [4] we can conclude that 
the probability of the event 

• both processes bave been coupled before time t or we have moved to stage b) before t 
is dominated by P(£i + £2 > tj. Continuing this line of reasoning we see that the coupling time 
Tcoupi is dominated by X!i=i & m the sense that 

P(r coupl >t)<p(53fi>i 



!=1 



Now we obtain 



j < e~ At E 









P(r COMp/ >t) <pfeCi 

= (E[e Afl ]) 5 e- At . 

Since E[e x ^] < 00 for all A < \{ a,b) , the claim follows. □ 

Remark 8. For pure probabilists the application of the more analytic Theorem^ via Corollary 
[^] may be somewhat unsatisfactorily; we included this result due to its applicability in the multi- 
dimensional setting. It seems to be very probable that Theorem^may be replaced by an additional 
coupling argument in the proof of Theorem 

6. Different jump distributions 

We end this work with some remarks concerning the case of a one-dimensional BM in (a, b) 
with two different jump distributions v a and v\,. It is a natural question, whether one can adopt 
the methods of this paper in order to prove another very recent results of Li, Leung and Rakesh 
[5] and Li and Leung [4], respectively, namely 

(6.1) sup 7i(i/„,i/(,) = Af 

v a ,v b 

and 

(6.2) inf 7i(i/ a , v b ) = Xq 

It is known that the supremum is (|6.ip is attained at v a = <5(( a +2&)/3 and Vb — ^(2a+b)/3 and that 
the infimum in (|6.2[) is never attained. Let us remark that equation (|6.2p is from a heuristic point 
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of view probabilistically rather clear, as taking v a;n := <5 a+ i/„ and i>b,n '■= ^b+i/n onc might expect 
that as n — > oo one gets a BM in (a, b) with reflecting boundary conditions. 

At the moment we can extend our methods only to certain very special classes of of jump 
distributions v a and v^, but these considerations indicate that with more effort one might be able 
to extend our methods to the more general situation of different jump distributions. Still at the 
moment we are far from being able to present a probabilistic proof of the beautiful assertions 
(|6.f|) and (|6.2p , but we believe that a solution of the following problem deepens the probabilistic 
understanding of the large time behavior of the BM with jump boundary considerably. 

Open Problem: Find a coupling approach to the recent results (|6.f [) and (|6.2p of Leung and Li. 

We hope to come back to this problem in a subsequent publication. 
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